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We reconstruct the radial profile of the projected gravitational potential of the galaxy cluster MACS J1206 from 592 spectroscopic 
measurements of velocities of cluster members. For doing so, we use a method we have developed recently based on the Richardson- 
Lucy deprojection algorithm and an inversion of the spherically-symmetric Jeans equation. We find that, within the uncertainties, our 
reconstruction agrees very well with a potential reconstruction from weak and strong gravitational lensing as well as with a potential 
obtained from X-ray measurements. In addition, our reconstruction is in good agreement with several common analytic profiles of 
the lensing potential. Varying the anisotropy parameter in the Jeans equation, we find that isotropy parameters which are either small, 
P < 0.2, or decrease with radius yield potential profiles which strongly disagree with that obtained from gravitational lensing. We 
achieve the best agreement between our potential profile and the profile from gravitational lensing if the anisotropy parameter rises 
quite steeply top x 0.6 within x: 0.5 Mpc and stays constant further out. 
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1. Introduction 

Galaxy clusters offer several classes of observables reflecting 
their overall internal constitution; gravitational lensing effects 
in their weak and strong variants. X-ray emission, the ther¬ 
mal Sunyaev-Zel’dovich (tSZ) effect and the kinematics of their 
member galaxies. We neglect radio emission, turbulence and the 
metal abundance in the intra cluster gas as well as the population 
statistics of the member galaxies here because these effects are 
locally driven. 

Gravitational lensing measures the gravitational tidal field, 
projected along the line-of-sight and thus directly probes the pro¬ 
jected gravitational potential, more precisely its curvature. The 
X-ray emission and the thermal Sunyaev-Zel’dovich effect de¬ 
pend on powers of the density and the temperature of the intra¬ 
cluster medium. If equilibrium assumptions hold, hydrostatic 
and virial equilibrium foremost, the X-ray emission, the thermal 
Sunyaev-Zel’dovich effect and galaxy kinematics will also be 
determined by the gravitational potential. Cmi'ent discussions of 
the validity of these equilibrium assumptions can be found e.g. 
Martino et al.| (|2014|l; von der Linden et al.| ( |2014[ ); |Planck| 


in 


ICollaboration et ah 1 2013 2014b|l. 


What is the gravitational potential that agrees best with all 
cluster observables? This question is relevant for different rea¬ 
sons. First, different observables trace the gravitational potential 
at different scales. The cluster core can be probed by the stel¬ 
lar kinematics of the brightest cluster galaxy, strong lensing and 
X-ray emission probe the innermost regions, weak lensing and 
galaxy kinematics probe large scales, and the tSZ effect falls in 
between (see Umetsu et al.| ( |2014l l; [Merten et al.|P015| l; [New-| 
|man et al.|(|2013|l). Aiming it a reliable reconstruction of cluster 


density profiles from their cores to their outskirts, combining all 
observables into a unique potential reconstruction offers the ad¬ 
vantage of covering all relevant scales in a single step. Second, 
comparing cluster potential reconstructions based on lensing on 
the one hand and based on the rest of the observables on the other 
hand allows testing the equilibrium assumptions or possible de¬ 
viations therefrom. Third, lensing and the other observables do 
not necessarily see the same gravitational potential. While lens¬ 
ing is sensitive to the sum of the Bardeen potentials, the other 
observables probe the spatial potential only. In general relativity, 
the two Bardeen potentials agree in case of negligible anisotropic 
stress. Differences in potential reconstructions based on lensing 
compared to other observables may also hint at deviations from 
relativity (see also Sartoris et al.||2014 Barreira et al.||20l3| in 
this context). 

Aiming at the gravitational potential has the major and im¬ 
portant advantages that it is a locally measureable quantity which 
is directly related to the observables listed above (see also | An- 1 
Igrick & Bartelmann|2009[[2012[ ). 

We have recently developed methods for reconstructing the 
projected gravitational potential of galaxy clusters from their X- 
ray emission, their tSZ effect and the kinematics of their member 
galaxies ( |Sarli et al.|2014[|Konrad et al.|2013[|Majer et al.|2013| l. 
They operate similarly, but with important differences in detail: 
An observable is deprojected by means of the Richardson-Lucy 
algorithm, requiring symmetry assumptions. The deprojected 
quantities are related to the three-dimensional gravitational po¬ 
tential by relations derived from justifiable equilibrium assump¬ 
tions. The gravitational potential can then be projected along 
the line-of-sight. These methods complement our techniques for 
joint cluster reconstruction from weak and strong gravitational 
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lensing (|Bartelmann et al. 1996 Cacciato et al. 12006[ 

Merten 

et al. 

20091 see also|Coe et al.|20121|Merten et al.|201 ![ 

Merten 

2014 

for examples). 



In this paper, we reconstruct the projected gravitational po¬ 
tential of the galaxy cluster MACS J1206.2-0847 based on 
galaxy kinematics, applying the technique developed in Sarli 
|et al.| ( |2OT4] l. This work is structured as follows: Sect. 2 briefly 
reviews the reconstruction method. In Sect. 3, we describe the 
selection and preparation of the data. Sect. 4 presents the re¬ 
sults and compares the projected gravitational potential to that 
obtained from gravitational lensing. We summarise in Sect. 5 
and discuss other methods in Sect. 6. 


2. Reconstruction method 


The basic assumptions of our reconstruction method are the fol¬ 
lowing: We treat the galaxy cluster as a spherically symmet¬ 
ric gas cloud of collisionless, pointlike test particles of mass 
m, i.e. galaxies, moving in the gravitational potential of their 
common dark matter halo. The system can then be described 
by the Jeans equation relating the radial velocity dispersion cr^ 
weighted by the galaxy number density pg^i to the gravitational 
potential (p (t>fm, 

1 g(Pgdcr^) , dcp 

Pgal Of f Of 


where the anisotropy parameter fi 1 - o-gfcr^ quantifies the 
ratio between the tangential and the radial velocity dispersions. 

Expanding on the formal analogy with gas dynamics, we de¬ 
fine an effective galaxy pressure P Pgaicr^. Using this def¬ 
inition and setting fi to zero, ([T]) would turn into the equation 
of hydrostatic equilibrium for a gas. In addition, we introduce a 
polytropic relation between the effective galaxy pressure and the 
matter density. 


P^Po 



( 2 ) 


This was justified for a simulated cluster in |Sarli et al.| ( |2014] ). In 
fact, the authors first tested with several density profiles that the 
polytropic assumption is reasonable. Additionally, we shall see 
in our reconstruction results that the exact choice of the poly¬ 
tropic index affects the final results only very mildly. 

The reconstruction algorithm described in detail in ||Sarli et al.| 
(|2014|l proceeds along the following steps (see also Fig.[]J: 


1. In order to obtain the effective pressure P, we first have to de¬ 
project the actual observable, i.e. the line-of-sight projected 
velocity dispersion weighted by the galaxy number density, 
PgaiCTj^s- This is done via the Richardson-Lucy deconvolution 
or deprojection, see [Lucy | ( | 1 974| [l994[ ). 

2. Furthermore, we make use of the polytropic relation •El to 
rewrite the Jeans equation Q in terms of the effective pres¬ 
sure. This leads to a Volterra integral equation of the sec¬ 
ond kind for the gravitational potential <t). After fixing the 
shape of the anisotropy profile, we can solve this equation in 
a quickly converging iteration process. 

3. Finally, we project <1) along the line-of-sight to find the pro¬ 
jected gravitational potential T'. 


This entire algorithm depends on three parameters and 
a function being carefully adapted during the reconstruction. 
Noise suppression within the Richardson-Lucy deprojection re¬ 
quires regularisation, controlled by a smoothing scale L and an 
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Fig. 1: Scheme of the reconstruction algorithm. The observable 
is converted to a three-dimensional quantity by Richardson-Lucy 
deprojection. This is then turned into the three-dimensional grav¬ 
itational potential, solving the Jeans equation. The resulting po¬ 
tential is Anally projected along the line-of-sight. 


amplitude a. Furthermore, we have to choose the polytropic in- 
dex y in ([ 2 1 , usually being of order unity, as cluster analyses in 
|Sarli et ah ( 2014| l suggest. This also implies that the galaxy fluid 
can approximately be treated as an isothermal gas. The impor¬ 
tance of the polytropic index is investigated in more detail in 
Sect. 4.1. 


As we shall demonstrate later, the dominant parameter is the 
anisotropy profile /l(r). However, if kinematic data are used ex¬ 
clusively for cluster reconstruction, there is in principle a de¬ 
generacy between the anisotropy parameter and the gravitational 
potential. We resolve this well-known degeneracy by fixing the 
yS-proflle in order to obtain a non-parametric estimate for the 
gravitational potential. To avoid an arbitrary, unmotivated guess 
for y8(r), we will choose it such that the reconstructed potential 
agrees best with a reconstruction based on gravitational-lensing 
data. Alternative methods for breaking the anisotropy-mass de¬ 
generacy and for recovering the gravitational potential using 
galaxy kinematics will be dicussed in the last section. 


3. Data selection and preparation 


The observational input data for the relaxed and massive galaxy 


cluster MACS J1206 (Miens ~ 5 ■ lO'^Mo/r-* [Urnetsu et al. 
P012| l) investigated here consist of a catalogue of member 
galaxies |Biviano et al.| (|2013|l, observed in the context of the 


CLASH project ( Postman et al.|[2012| l as part of a large spec 
troscopic campaign carried out with the Very Large Telescope 
(VLT) (CLASH-VLT Large Programme; |Rosati et al. p014|l). 


According to Biviano et al. (2013|l and Annunziate’ 


la et al. 


P014| l, the spatial incompleteness of the spectroscopic sample 
varies with position in the cluster by less than 20%. 

The following analysis is based on 592 member galaxies in 
total, selected via the method presented in |Biviano et al.| ( |2013| l. 
In order to arrive at the line-of-sight projected velocity disper¬ 
sion profile, we assign to every cluster member a projected ra¬ 
dius from the cluster centre, which we take to be marked by the 
brightest cluster galaxy. This is achieved by multiplying the an¬ 
gular separation with the angular-diameter distance, calculated 
with a standard ACDM-cosmology with recent cosmological pa¬ 
rameters from |Planck Collaboration et al.| ( |2014a| l. 

We determine the line-of-sight projected velocities from the 
measured redshifts as described in |Harrison| ( |1974| l. In our case 
we can use the non-relativistic relation between redshift z and 
line-of-sight velocity yios, since all involved velocities (see be¬ 
low) are clearly in the non-relativistic regime. Thus, 


^los 


cz . 


(3) 
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Due to the gauge freedom of the lensing potential, one is in 
particular allowed to add an arbitrary constant. By convention, 
we scale the potential such that 'l'(O) = 1 and use the normalisa¬ 
tion 'I'(/?cut) = 0 for a given cut-off radius Rcut = 3 Mpc. This is 
even larger than the virial radius of MACS J1206, being approx¬ 
imately 2 Mpc ( |Umetsu et al.|2012] l. 

In order to obtain error bars and to suppress the effect of 
outliers, we perform a bootstrap analysis, i.e. for each bootstrap 
sample, we draw as many times with replacement from the orig¬ 
inal dataset as it has member galaxies. Then, we reconstruct the 
gravitational potential for each individual bootstrap sample. This 
procedure is repeated 300 times. Finally, we calculate the mean 
of all potentials and their standard deviation. The bootstrap anal¬ 
ysis also allows to assess the effects of incompleteness and un¬ 
certain membership assignment by varying the sample of galax¬ 
ies entering into the potential reconstruction. 


Fig. 2: Velocities projected along the l.o.s. of cluster members. 


As we are only interested in the galaxy velocities with respect 
to the cluster’s centre-of-mass, we have to take into account that 
the observed redshift includes four main contributions: the mo¬ 
tion of the observer with respect to a local comoving observer, 
the motion of the cluster’s centre-of-mass with respect to a local 
comoving observer, the motion of each galaxy with respect to the 


cluster’s centre-of-mass, and finally the Hubble expansion. Har- 
|rison| ( [T974) l has shown that neglecting the motion of the cluster’s 
centre-of-mass leads to a multiplication of redshifts. 


1 + Z - (1 + 2o)(l + 2cosm)(l + Zg) > 


(4) 


where zo is the redshift due to the observer’s motion with respect 
to a local comoving observer, Zcosm represents the Hubble ex¬ 
pansion and zg the motion of each galaxy with respect to a local 
comoving observer. To leading order, we can calculate the clus¬ 
ter average of Q and use that the average redshift of all member 
galaxies with respect to the local comoving frame vanishes be¬ 
cause we assume an isotropic distribution of velocities along the 
l.o.s.. 


(ZG> = 0 . 


(5) 


Furthermore, we can safely assume that the contribution by the 
Hubble expansion is the same for all cluster members. 


{^cosm) — ^cosm ■ 

Thus, we end up with 

l+<z> 


1 T ^cosm — 


1 +Z0 


(6) 


(7) 


Inserting 0 into we find the following expression for wios, 
see also figure]^ 


0OS “ ^ Zg ~ ^ 


1 + z 

1+<Z) 


-1 


( 8 ) 


To arrive at the density-weighted, projected velocity disper¬ 
sion profile Pgaicrfos, we adapt the bin width to the one chosen for 
the potential from lensing for better comparability, i.e. a constant 
bin width of 0.1 Mpc. Next, we calculate the velocity dispersion 
within each bin and reconstruct the projected gravitational po¬ 
tential as described above by solving the radial Jeans equation. 


4. Results 

4. 1. Effect of reconstruction parameters 

As indicated above, we have to fix three reconstruction param¬ 
eters in order to perform the reconstruction: the regularisation 
amplitude a, the smoothing scale L and the polytropic index y. 
Additionally, we have to model the anisotropy profile /J(r). 

Figure[^shows how the variation of these reconstruction pa¬ 
rameters affects the resulting potential T*. The amplitude of the 
anisotropy parameter yS clearly has the dominant effect. In par¬ 
ticular, an isotropic velocity distribution corresponding to ^ 0 
seems very implausible because of the strongly varying shape of 
the potential resulting from this assumption. 

Although the assumption of a poly tropic relation Q may ap¬ 
pear quite bold. Fig. shows that the particular choice of the 
polytropic index does not really matter. Also the effects of the 
regularisation parameters a and L remain within the uncertain¬ 
ties, being of similar order as shown in Fig. When they are 
not varied, the parameters are set as follows: a - 0.3, L - 0.6, 
7=1.1 and [i - 0.6. The value chosen for the polytropic index y 
is motivated by the result of |Sarli et al.|p014j l that y would usu¬ 
ally be around unity. The regularisation parameters are chosen 
such that the agreement with the lensing reconstruction is best. 

4.2. Variable anisotropy profile 

We also systematically investigate the effect of an anisotropy pa¬ 
rameter varying with radius. For doing so, we compare radial 
profiles of gravitational potentials obtained with an anisotropy 
parameter increasing or decreasing linearly in discrete steps ev¬ 
ery 0.5 Mpc. In case of the decreasing anisotropy-parameter pro¬ 
file, we choose p - 0.8, 0.7, 0.6, 0.5, 0.4, 0.3, and reverse for the 
increasing case. 

Figure]^ shows both cases together with the potential recon¬ 
struction from weak and strong gravitational lensing and a re¬ 
construction with constant p - 0.6. First, one can observe that 
the qualitative form of the anisotropy profile affects the curvature 
of the potential: An increasing p profile leads to a potential with 
negative curvature. In contrast, a decreasing /J-profile leads to a 
positive curvature and strong variations in the potential at radii 
beyond 1 Mpc. The reconstruction using a constant p leads to an 
almost straight line for the radial potential profile. 

These results motivate the shape of the anisotropy profile 
used later in 0. Since the potential from lensing first shows 
a negative curvature at small radii and a linear evolution later. 
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Fig. 3; Reconstructed radial profiles of the projected gravitational potential of MACS J1206, obtained using different parameter sets 
for the smoothing amplitude a (panel a), smoothing scale L (panel b), the polytropic index y (panel c) and the anisotropy parameter 
j3 (panel d). 


we choose an anisotropy profile increasing for small radii and 
turning constant afterwards to reproduce the linear behaviour. 


4.3. Comparison with iensing and X-ray data 


As a first test, we compare the projected gravitational potential 
reconstructed from galaxy velocities with a reconstruction based 
on weak and strong gravitational Iensing data of [Merten et al.| 
( 2015| l using constraints from Zitrin et al. (20121. After convert¬ 
ing the 2D-lensing potential map (see figure]^ into a radial pro¬ 
file and using the same normalisation process as above, we can 
compare the result from gravitational Iensing to the potential ob¬ 
tained from cluster kinematics (Fig.|^. 


The reconstruction parameters are chosen as listed above. 
For the radial profile of the anisotropy parameter j3, we assume 
the following radially dependent profile, being motivated in Sect. 
4.2. This profile turns out to be in qualitative agreement with the 
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one in|Biviano et al.|(|2013|l: 


Pir) = 


0.1 

0.4 


0.6 


r <0.2 

0.2 < r < 0.5 . 
0.5 <r< 3.0 


(9) 


Additionally, we compare it with a reconstruction based on X- 
ray emission of the intracluster gas based on [Donahue et al.| 
P014] l using Chandra data up to 0.72 Mpc, which is a typical 
range of validity using this method for galaxy clusters. 

The growing difference between the Iensing and kinematic pro¬ 
file outside ~ 2 Mpc can be explained by the limited range of 
validity of each reconstruction method. The small difference at 
the innermost radii is due to increasing baryonic effects causing 
the equilibrium assumption ultimately to break down. The X-ray 
data agree quite well within the error boundaries. 


















































Dennis Stock et al.: The projected gravitational potential of the galaxy cluster MACS J1206 derived from galaxy kinematics 


CO 



Projected radius s [Mpc] 

Fig. 4; Radial profiles of the projected gravitational potential, 
reconstructed from galaxy kinematics with different prohles of 
the anisotropy parameter fi. The other reconstruction parameters 
remain hxed at a = 0.3, L - 0.6 and y - 1.1. 



Fig. 5: Two dimensional map of the lensing potential of 
MACS 1206 based on weak and strong lensing. 


4.4. Comparison with common iensing potentiais 

Furthermore, we can compare the radial prohle of the projected 
gravitational potential obtained from galaxy kinematics with ra¬ 
dial prohles common in gravitational-lensing studies, such as the 
singular (SIS) and the non-singular (NIS) isothermal spheres and 
the NFW prohle ([Navarro et al.||1996[ jNarayan & Bartelmannj 
19^ ; see Fig.0 

'Tsis(s) = A ■ s + B 

'Tnis(s) = C ■ ^sl + s^ + D 

'Pnfw('S) = E ■ h(x) + F (10) 


with the function h(x) depending on the radius in units of the 
scale radius .r := rjr^ (seejGolse & Kneib|p002|l): 


h{x) - < 


X 1 

In^ - - arcosh^- (x < 1) 

2 X 

X 1 

In^ - -H arccos^ - (x > 1) 

2 X 


( 11 ) 


where Sc, rs and A . ..F are free parameters being determined 
by htting the models ( [T0) i to the potential obtained from galaxy 
kinematics. The parameter values are summarised in Tabj^. 

We can compare the results for the scales Sc in case of the 
softened isothermal sphere and rj in case of NFW with |Umetsu| 
|et al.| ( |2OT2] l. They obtain: 

Sc = 0.000028 + 0.000006 (12) 

rs = 0.345 + 0.050 . (13) 


The agreement is quantihed by a goodness-of-ht parameter Q, 
which we dehne in analogy to the -function appropriate for 
uncorrelated measurements. Since our data points are not inde¬ 
pendent but correlated by the bootstrap method, i.e. some data 
points appear multiple times, we cannot interpret Q as a 
function. Because we just want to single out the best model, we 
do not go into a more elaborate analysis for correlated measure¬ 
ments here. For the /-th data point x, with standard deviation cr, 
and the corresponding model prediction /(x,), we have 



Fig. 6: Radial prohles of the gravitational potential of 
MACS J1206, reconstructed from combined strong and weak 
gravitational lensing (blue). X-ray emission (yellow) and from 
galaxy kinematics (red) with the parameters a = 0.3, L - 0.6, 
y - andjS as in (|^. 



As one can already see by eye, the softened isothermal sphere 
and the NFW-model describe the data equally well. More quan¬ 
titatively, the goodness-of-ht parameters for all lensing model 
are: 

^sis = 6.09 
f^Nis = 4.00 

f^NFw — 3.54 . (15) 

So far, all analytic lensing potentials are compatible with the re¬ 
construction from galaxy kinematics, but the error bars are sub¬ 
stantial. 

5. Summary 

We have recently developed a method for reconstructing the pro¬ 
jected gravitational potential of galaxy clusters from the kine¬ 
matics of their member galaxies. In the study described in this 
paper, we have applied this method for the hrst time to a real 
galaxy cluster (MACS J1206). The reconstruction assumes that 
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Table 1: Best fit parameters for different models 


Model 

Parameters 

Singular Isothermal Sphere 

A = -0.35 + 0.01 

B = 1.02 + 0.01 


Non-singular Isothermal Sphere 

C = -0.39 + 0.01 

Sc = 0.38 + 0.06 

D = 1.11+0.02 

NEW 

E = -0.84 + 0.10 

r, = 1.37 + 0.15 

F = 0.97 + 0.01 



Fig. 7; Comparison of three analytic profiles of gravitational- 
lensing potentials with the gravitational-potential profile of 
MACS J1206 reconstructed from galaxy kinematics using a - 
0.3, L - 0.6, 7=1.1 and /? as in (|^. 


the three-dimensional gravitational potential of this cluster is ap¬ 
proximately spherically symmetric. The potential reconstruction 
is based on 592 measured galaxy velocities. 

Our results can be summarised as follows: 

- The reconstruction algorithm requires four parameters to be 
set: the anisotropy parameter /?; the polytropic index y of 
the effective galaxy pressure; the amplitude a of the regular- 
isation term in the Richardson-Lucy deconvolution and the 
length scale L of the regularisation. We found that the recon¬ 
structed potential is almost insensitive to the “nuisance pa¬ 
rameters” 7 , a and L, while it does depend quite sensitively 
on the anisotropy parameter /?. 

- Within the bootstrap error bars, the reconstructed potential 
profile is indistinguishable from the potential profile ob¬ 
tained from the combination of weak and strong gravita¬ 
tional lensing and from X-ray analysis. Typical analytic lens- 
ing profiles, such as the singular and non-singular isothermal 
spheres and the NFW profile also agree well with our recon¬ 
struction from galaxy kinematics. 

- Small anisotropy parameters y6 < 0.2 or anisotropy- 

parameter profiles decreasing with radius yield potential pro¬ 
files differing strongly from the profile obtained from gravi¬ 
tational lensing. Best agreement with gravitational lensing is 
achieved with an anisotropy parameter increasing to y6 ^ 0.6 
within 0.5 Mpc and staying constant further out. 

Our results show that the method returns convincing results 
on one well-studied galaxy cluster which we use here as a test 
case. Even though the uncertainties are still large, the anisotropy 
parameter can be rather well constrained by comparing our re¬ 
sults with those obtained from gravitational lensing. We see three 
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main future applications of this method: (1) in joint reconstruc¬ 
tions of galaxy-cluster potentials compatible with all cluster ob¬ 
servables; (2) in constraints of the anisotropy parameter; and (3) 
in tests of fundamental assumptions such as hydrostatic or virial 
equilibrium, and possibly also of general relativity, if applied to 
large and well-measured cluster samples. 


6. Discussion of other methods 


On top of the conditions leading to the Jeans equation, the 
method presented above implicitly rests on the assumption that 
all test particles (i.e. the galaxies) have the same mass, that their 
number density is proportional to their mass density and that one 
can establish a polytropic relation between their mass density p 
and the effective galaxy pressure P := Pgaicr^. Fixing the ve¬ 
locity anisotropy profile y6(r) enables us to constrain the gravita¬ 
tional potential of the cluster via the radial Jeans equation using 
the observed density-weighted galaxy velocity dispersion. Our 
method does not assume the mass to follow light a priori, i.e. 
Pdm does not necessarily need to follow pg^i. Since the resulting 
potential is non-parametric, one can test the validity of particular 
parametric models as done above. 

Another common deprojection method is the Abel inversion 
using derivatives of observables. However, due to the fluctuat¬ 
ing nature of most observables which cause strongly fluctuating 
derivatives, we decided to use the Richardson-Lucy deprojection 
instead which only involves integrals. 

There are many other methods to determine the gravitational 
potential of a cluster, or its mass, or to constrain the velocity 
anisotropy. For example, Binney & Mamon ( |1982| l use the line- 
of-sight velocity dispersions as observables together with the 
surface brightness in spherical galaxies. Assuming a constant 
mass-to-light ratio, they can replace the unknown density in the 
Jeans equation by the luminosity, and can thus determine the ve¬ 
locity anisotropy profile /3{r). However, in the context of galaxy 
clusters it is not clear that the same relation between mass and 
luminosity should hold as well. 

[Dejonghe & Merri'tt| ( |1992) l determine the potential <1) by tak¬ 
ing into account higher than second-order velocity moments of 
the collision-less Boltzmann equation (CBE). Given the velocity 
dispersions, and making use of a linear relation between the ob¬ 
servables and O, they expand the gravitational potential in terms 
of basis functions with unspecified coefficients. Using the posi¬ 
tivity constraint of the distribution function and higher order mo¬ 
ments of the CBE, they are able to determine the coefficients and 
thus the potential. However, from a practical point of view, it be¬ 
comes increasingly difficult to constrain higher-order moments 
of the velocity distribution from their observable line-of-sight 
projections, as the projection integrals turn out to be substan¬ 
tially more complex with increasing order (see [Richardson &| 
|Fairbairn](|2013|l for instance). 


A method very similar to ours is proposed by Mamon & 
|Boue| ( [2010| l: starting again from the Jeans equation and fixing 
yS(r), they determine both the number density and the velocity 
dispersion observationally. Besides their using the Abel inver¬ 
sion, the only difference to our method is that they determine 
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the number density of galaxies appearing in the Jeans equation 
by counting. Hence, they do not need our assumed polytropic 
relation between density and effective pressure. In future work, 
one could use this to test directly to what degree the polytropic 
assumption is appropriate. 

[Wolf et ar] ( |2010[ ) and |van der Marel et al.| ( |2000l l are consid¬ 
ering a general parametrisation of /3(r) and of the mass M(r) and 
constrain all free parameters via a maximum-likelihood analy¬ 
sis of the velocity dispersion predicted by the Jeans equation 
compared to the observed one, [van der Marel et al.|(2000| l use 
a constant velocity anisotropy whereas Wolf et al. ( |2()10[ ) take a 
varying parametric profile. 

Another way to break the anisotropy mass degeneracy is 
taking the fourth order velocity moments of the CBE into ac¬ 
count ( |Lokas| ( |2002| l, see also [Richardson & Fairbairn] ( |2013| l in 
this context). By expressing the velocity distribution function in 
terms of two integrals of motion, the energy E and the angular 
momentum L, and further assuming that the distribution is sepa¬ 
rable in E and L with fi - const., the two fourth-order moments 
are reduced to one equation involving the radial velocity to the 
fourth power, (3, the radial velocity dispersion and <D>. Thus, by 
measuring the velocity dispersion and one can infer the grav¬ 
itational potential by solving this fourth-order equation together 
with the Jeans equation. 

The MAMPOSSt method introduced by jMamon et al.|p013| l 
breaks the degeneracy assuming parametric forms for the grav¬ 
itational potential, the velocity anisotropy and the distribution 
of 3D-velocities. By performing a maximum-likelihood fit of 
the galaxy distribution in the projected phase space, all param¬ 
eters can be determined. However, it has recently been shown 
( [Richardson & Fairbairn| ( |2013| )) that their assumption of a Gaus¬ 
sian 3D-velocity distribution is incompatible with the equilib¬ 
rium assumption underlying the CBF. One would thus have to 
choose a different 3D-velocity distribution. 
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